Spaceflight effects on human vascular smooth muscle cell phenotype and function.

The cardiovascular system is strongly impacted by the hazards of spaceflight. Astronauts spending steadily increasing lengths of time in microgravity are subject to cardiovascular deconditioning resulting in loss of vascular tone, reduced total blood volume, and diminished cardiac output. Appreciating the mechanisms by which the cells of the vasculature are altered during spaceflight will be integral to understanding and combating these deleterious effects as the human presence in space advances. In this study, we performed RNA-Seq analysis coupled with review by QIAGEN Ingenuity Pathway Analysis software on human aortic smooth muscle cells (HASMCs) cultured for 3 days in microgravity and aboard the International Space Station to assess the transcriptomic changes that occur during spaceflight. The results of our RNA-Seq analysis show that SMCs undergo a wide range of transcriptional alteration while in space, significantly affecting 4422 genes. SMCs largely down-regulate markers of the contractile, synthetic, and osteogenic phenotypes including smooth muscle alpha actin (αSMA), matrix metalloproteinases (MMPs), and bone morphogenic proteins (BMPs). Additionally, components of several cellular signaling pathways were strongly impacted including the STAT3, NFκB, PI3K/AKT, HIF1α, and Endothelin pathways. This study highlights the significant changes in transcriptional behavior SMCs exhibit during spaceflight and puts these changes in context to better understand vascular function in space.

N6-methyladenosine in 7SK small nuclear RNA underlies RNA polymerase II transcription regulation.

N6-methyladenosine (m6A) modifications play crucial roles in RNA metabolism. How m6A regulates RNA polymerase II (RNA Pol II) transcription remains unclear. We find that 7SK small nuclear RNA (snRNA), a regulator of RNA Pol II promoter-proximal pausing, is highly m6A-modified in non-small cell lung cancer (NSCLC) cells. In A549 cells, we identified eight m6A sites on 7SK and discovered methyltransferase-like 3 (METTL3) and alkB homolog 5 (ALKBH5) as the responsible writer and eraser. When the m6A-7SK is specifically erased by a dCasRx-ALKBH5 fusion protein, A549 cell growth is attenuated due to reduction of RNA Pol II transcription. Mechanistically, removal of m6A leads to 7SK structural rearrangements that facilitate sequestration of the positive transcription elongation factor b (P-TEFb) complex, which results in reduction of serine 2 phosphorylation (Ser2P) in the RNA Pol II C-terminal domain and accumulation of RNA Pol II in the promoter-proximal region. Taken together, we uncover that m6A modifications of a non-coding RNA regulate RNA Pol II transcription and NSCLC tumorigenesis.

Myotonic dystrophy type 1 embryonic stem cells show decreased myogenic potential, increased CpG methylation at the DMPK locus and RNA mis-splicing.

Skeletal muscle tissue is severely affected in myotonic dystrophy type 1 (DM1) patients, characterised by muscle weakness, myotonia and muscle immaturity in the most severe congenital form of the disease. Previously, it was not known at what stage during myogenesis the DM1 phenotype appears. In this study we differentiated healthy and DM1 human embryonic stem cells to myoblasts and myotubes and compared their differentiation potential using a comprehensive multi-omics approach. We found myogenesis in DM1 cells to be abnormal with altered myotube generation compared to healthy cells. We did not find differentially expressed genes between DM1 and non-DM1 cell lines within the same developmental stage. However, during differentiation we observed an aberrant inflammatory response and increased CpG methylation upstream of the CTG repeat at the myoblast level and RNA mis-splicing at the myotube stage. We show that early myogenesis modelled in hESC reiterates the early developmental manifestation of DM1.

Generation of a Novel SARS-CoV-2 Sub-genomic RNA Due to the R203K/G204R Variant in Nucleocapsid: Homologous Recombination has Potential to Change SARS-CoV-2 at Both Protein and RNA Level.

BACKGROUND: Genetic variations across the SARS-CoV-2 genome may influence transmissibility of the virus and the host's anti-viral immune response, in turn affecting the frequency of variants over time. In this study, we examined the adjacent amino acid polymorphisms in the nucleocapsid (R203K/G204R) of SARS-CoV-2 that arose on the background of the spike D614G change and describe how strains harboring these changes became dominant circulating strains globally. METHODS: Deep-sequencing data of SARS-CoV-2 from public databases and from clinical samples were analyzed to identify and map genetic variants and sub-genomic RNA transcripts across the genome. Results: Sequence analysis suggests that the 3 adjacent nucleotide changes that result in the K203/R204 variant have arisen by homologous recombination from the core sequence of the leader transcription-regulating sequence (TRS) rather than by stepwise mutation. The resulting sequence changes generate a novel sub-genomic RNA transcript for the C-terminal dimerization domain of nucleocapsid. Deep-sequencing data from 981 clinical samples confirmed the presence of the novel TRS-CS-dimerization domain RNA in individuals with the K203/R204 variant. Quantification of sub-genomic RNA indicates that viruses with the K203/R204 variant may also have increased expression of sub-genomic RNA from other open reading frames. CONCLUSIONS: The finding that homologous recombination from the TRS may have occurred since the introduction of SARS-CoV-2 in humans, resulting in both coding changes and novel sub-genomic RNA transcripts, suggests this as a mechanism for diversification and adaptation within its new host.

Generation of a novel SARS-CoV-2 sub-genomic RNA due to the R203K/G204R variant in nucleocapsid: homologous recombination has potential to change SARS-CoV-2 at both protein and RNA level.

BACKGROUND: Genetic variations across the SARS-CoV-2 genome may influence transmissibility of the virus and the host’s anti-viral immune response, in turn affecting the frequency of variants over-time. In this study, we examined the adjacent amino acid polymorphisms in the nucleocapsid (R203K/G204R) of SARS-CoV-2 that arose on the background of the spike D614G change and describe how strains harboring these changes became dominant circulating strains globally. METHODS: Deep sequencing data of SARS-CoV-2 from public databases and from clinical samples were analyzed to identify and map genetic variants and sub-genomic RNA transcripts across the genome. RESULTS: Sequence analysis suggests that the three adjacent nucleotide changes that result in the K203/R204 variant have arisen by homologous recombination from the core sequence (CS) of the leader transcription-regulating sequence (TRS) rather than by stepwise mutation. The resulting sequence changes generate a novel sub-genomic RNA transcript for the C-terminal dimerization domain of nucleocapsid. Deep sequencing data from 981 clinical samples confirmed the presence of the novel TRS-CS-dimerization domain RNA in individuals with the K203/R204 variant. Quantification of sub-genomic RNA indicates that viruses with the K203/R204 variant may also have increased expression of sub-genomic RNA from other open reading frames. CONCLUSIONS: The finding that homologous recombination from the TRS may have occurred since the introduction of SARS-CoV-2 in humans resulting in both coding changes and novel sub-genomic RNA transcripts suggests this as a mechanism for diversification and adaptation within its new host.

RNA structure probing to characterize RNA-protein interactions on a low abundance pre-mRNA in living cells.

In vivo RNA structure analysis has become a powerful tool in molecular biology, largely due to the coupling of an increasingly diverse set of chemical approaches with high-throughput sequencing. This has resulted in a transition from single target to transcriptome-wide approaches. However, these methods require sequencing depths that preclude studying low abundance targets, which are not sufficiently captured in transcriptome-wide approaches. Here we present a ligation-free method to enrich for low abundance RNA sequences, which improves the diversity of molecules analyzed and results in improved analysis. In addition, this method is compatible with any choice of chemical adduct or read-out approach. We utilized this approach to study an autoregulated event in the pre-mRNA of the splicing factor, muscleblind-like splicing regulator 1 (MBNL1).

Cell-type-specific dysregulation of RNA alternative splicing in short tandem repeat mouse knockin models of myotonic dystrophy.

Short tandem repeats (STRs) are prone to expansion mutations that cause multiple hereditary neurological and neuromuscular diseases. To study pathomechanisms using mouse models that recapitulate the tissue specificity and developmental timing of an STR expansion gene, we used rolling circle amplification and CRISPR/Cas9-mediated genome editing to generate Dmpk CTG expansion (CTGexp) knockin models of myotonic dystrophy type 1 (DM1). We demonstrate that skeletal muscle myoblasts and brain choroid plexus epithelial cells are particularly susceptible to Dmpk CTGexp mutations and RNA missplicing. Our results implicate dysregulation of muscle regeneration and cerebrospinal fluid homeostasis as early pathogenic events in DM1.

RAN Translation Regulated by Muscleblind Proteins in Myotonic Dystrophy Type 2.

Several microsatellite-expansion diseases are characterized by the accumulation of RNA foci and RAN proteins, raising the possibility of a mechanistic connection. We explored this question using myotonic dystrophy type 2, a multisystemic disease thought to be primarily caused by RNA gain-of-function effects. We demonstrate that the DM2 CCTG⋅CAGG expansion expresses sense and antisense tetrapeptide poly-(LPAC) and poly-(QAGR) RAN proteins, respectively. In DM2 autopsy brains, LPAC is found in neurons, astrocytes, and glia in gray matter, and antisense QAGR proteins accumulate within white matter. LPAC and QAGR proteins are toxic to cells independent of RNA gain of function. RNA foci and nuclear sequestration of CCUG transcripts by MBNL1 is inversely correlated with LPAC expression. These data suggest a model that involves nuclear retention of expansion RNAs by RNA-binding proteins (RBPs) and an acute phase in which expansion RNAs exceed RBP sequestration capacity, are exported to the cytoplasm, and undergo RAN translation. VIDEO ABSTRACT.

Characterization of the UGA-recoding and SECIS-binding activities of SECIS-binding protein 2.

Selenium, a micronutrient, is primarily incorporated into human physiology as selenocysteine (Sec). The 25 Sec-containing proteins in humans are known as selenoproteins. Their synthesis depends on the translational recoding of the UGA stop codon to allow Sec insertion. This requires a stem-loop structure in the 3' untranslated region of eukaryotic mRNAs known as the Selenocysteine Insertion Sequence (SECIS). The SECIS is recognized by SECIS-binding protein 2 (SBP2) and this RNA:protein interaction is essential for UGA recoding to occur. Genetic mutations cause SBP2 deficiency in humans, resulting in a broad set of symptoms due to differential effects on individual selenoproteins. Progress on understanding the different phenotypes requires developing robust tools to investigate SBP2 structure and function. In this study we demonstrate that SBP2 protein produced by in vitro translation discriminates among SECIS elements in a competitive UGA recoding assay and has a much higher specific activity than bacterially expressed protein. We also show that a purified recombinant protein encompassing amino acids 517-777 of SBP2 binds to SECIS elements with high affinity and selectivity. The affinity of the SBP2:SECIS interaction correlated with the ability of a SECIS to compete for UGA recoding activity in vitro. The identification of a 250 amino acid sequence that mediates specific, selective SECIS-binding will facilitate future structural studies of the SBP2:SECIS complex. Finally, we identify an evolutionarily conserved core cysteine signature in SBP2 sequences from the vertebrate lineage. Mutation of multiple, but not single, cysteines impaired SECIS-binding but did not affect protein localization in cells.

Alternative transcripts and 3'UTR elements govern the incorporation of selenocysteine into selenoprotein S.

Selenoprotein S (SelS) is a 189 amino acid trans-membrane protein that plays an important yet undefined role in the unfolded protein response. It has been proposed that SelS may function as a reductase, with the penultimate selenocysteine (Sec(188)) residue participating in a selenosulfide bond with cysteine (Cys(174)). Cotranslational incorporation of Sec into SelS depends on the recoding of the UGA codon, which requires a Selenocysteine Insertion Sequence (SECIS) element in the 3'UTR of the transcript. Here we identify multiple mechanisms that regulate the expression of SelS. The human SelS gene encodes two transcripts (variants 1 and 2), which differ in their 3'UTR sequences due to an alternative splicing event that removes the SECIS element from the variant 1 transcript. Both transcripts are widely expressed in human cell lines, with the SECIS-containing variant 2 mRNA being more abundant. In vitro experiments demonstrate that the variant 1 3'UTR does not allow readthrough of the UGA/Sec codon. Thus, this transcript would produce a truncated protein that does not contain Sec and cannot make the selenosulfide bond. While the variant 2 3'UTR does support Sec insertion, its activity is weak. Bioinformatic analysis revealed two highly conserved stem-loop structures, one in the proximal part of the variant 2 3'UTR and the other immediately downstream of the SECIS element. The proximal stem-loop promotes Sec insertion in the native context but not when positioned far from the UGA/Sec codon in a heterologous mRNA. In contrast, the 140 nucleotides downstream of the SECIS element inhibit Sec insertion. We also show that endogenous SelS is enriched at perinuclear speckles, in addition to its known localization in the endoplasmic reticulum. Our results suggest the expression of endogenous SelS is more complex than previously appreciated, which has implications for past and future studies on the function of this protein.

Identification of a signature motif for the eIF4a3-SECIS interaction.

eIF4a3, a DEAD-box protein family member, is a component of the exon junction complex which assembles on spliced mRNAs. The protein also acts as a transcript-selective translational repressor of selenoprotein synthesis during selenium deficiency. Selenocysteine (Sec) incorporation into selenoproteins requires a Sec Insertion Sequence (SECIS) element in the 3' untranslated region. During selenium deficiency, eIF4a3 binds SECIS elements from non-essential selenoproteins, preventing Sec insertion. We identified a molecular signature for the eIF4a3-SECIS interaction using RNA gel shifts, surface plasmon resonance and enzymatic foot printing. Our results support a two-site interaction model, where eIF4a3 binds the internal and apical loops of the SECIS. Additionally, the stability of the complex requires uridine in the SECIS core. In terms of protein requirements, the two globular domains of eIF4a3, which are connected by a linker, are both critical for SECIS binding. Compared to full-length eIF4a3, the two domains in trans bind with a lower association rate but notably, the uridine is no longer important for complex stability. These results provide insight into how eIF4a3 discriminates among SECIS elements and represses translation.

Eukaryotic initiation factor 4a3 is a selenium-regulated RNA-binding protein that selectively inhibits selenocysteine incorporation.

The synthesis of selenoproteins requires the translational recoding of the UGA stop codon as selenocysteine. During selenium deficiency, there is a hierarchy of selenoprotein expression, with certain selenoproteins synthesized at the expense of others. The mechanism by which the limiting selenocysteine incorporation machinery is preferentially utilized to maintain the expression of essential selenoproteins has not been elucidated. Here we demonstrate that eukaryotic initiation factor 4a3 (eIF4a3) is involved in the translational control of a subset of selenoproteins. The interaction of eIF4a3 with the selenoprotein mRNA prevents the binding of SECIS binding protein 2, which is required for selenocysteine insertion, thereby inhibiting the synthesis of the selenoprotein. Furthermore, the expression of eIF4a3 is regulated in response to selenium. Based on knockdown and overexpression studies, eIF4a3 is necessary and sufficient to mediate selective translational repression in cells. Our results support a model in which eIF4a3 links selenium status with differential selenoprotein expression.

Known turnover and translation regulatory RNA-binding proteins interact with the 3' UTR of SECIS-binding protein 2.

The human selenoproteome is composed of approximately 25 selenoproteins, which cotranslationally incorporate selenocysteine, the 21st amino acid. Selenoprotein expression requires an unusual translation mechanism, as selenocysteine is encoded by the UGA stop codon. SECIS-binding protein 2 (SBP2) is an essential component of the selenocysteine insertion machinery. SBP2 is also the only factor known to differentiate among selenoprotein mRNAs, thereby modulating the relative expression of the individual selenoproteins. Here, we show that expression of SBP2 protein varies widely across tissues and cell types examined, despite previous observations of only modest variation in SBP2 mRNA levels. This discrepancy between SBP2 mRNA and protein levels implies translational regulation, which is often mediated via untranslated regions (UTRs) in regulated transcripts. We have identified multiple sequences in the SBP2 3' UTR that are highly conserved. The proximal short conserved region is GU rich and was subsequently shown to be a binding site for CUG-BP1. The distal half of the 3' UTR is largely conserved, and multiple proteins interact with this region. One of these proteins was identified as HuR. Both CUG-BP1 and HuR are members of the Turnover and Translation Regulatory RNA-Binding Protein family (TTR-RBP). Members of this protein family are linked by the common ability to rapidly effect gene expression through alterations in the stability and translatability of target mRNAs. The identification of CUG-BP1 and HuR as factors that bind to the SBP2 3' UTR suggests that TTR-RBPs play a role in the regulation of SBP2, which then dictates the expression of the selenoproteome.

Altered RNA binding activity underlies abnormal thyroid hormone metabolism linked to a mutation in selenocysteine insertion sequence-binding protein 2.

The expression of selenoproteins requires the translational recoding of the UGA stop codon to selenocysteine. In eukaryotes, this requires an RNA stem loop structure in the 3'-untranslated region, termed a selenocysteine insertion sequence (SECIS), and SECIS-binding protein 2 (SBP2). This study implicates SBP2 in dictating the hierarchy of selenoprotein expression, because it is the first to show that SBP2 distinguishes between SECIS elements in vitro. Using RNA electrophoretic mobility shift assays, we demonstrate that a naturally occurring mutation in SBP2, which correlates with abnormal thyroid hormone function in humans, lies within a novel, bipartite RNA-binding domain. This mutation alters the RNA binding affinity of SBP2 such that it no longer stably interacts with a subset of SECIS elements. Assays performed under competitive conditions to mimic intracellular conditions suggest that the differential affinity of SBP2 for various SECIS elements will determine the expression pattern of the selenoproteome. We hypothesize that the selective loss of a subset of selenoproteins, including some involved in thyroid hormone homeostasis, is responsible for the abnormal thyroid hormone metabolism previously observed in the affected individuals.